Josephson effect in superfiuid atomic Fermi-gases 
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We consider an analog of the internal Josephson effect in superfiuid atomic Fermi-gases. Four 
different hyperfine states of the atoms are assumed to be trapped and to form two superfluids 
via the BCS-type of pairing. We show that Josephson oscillations can be realized by coupling 
the superfluids with two laser fields. Choosing the laser detunings in a suitable way leads to an 
asymmetric below-gap tunneling effect for which there exists no analogue in the context of solid-state 
superconductivity. 
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Cooling of trapped gases of Fermionic atoms well be- 
low the Fermi temperature [j], ||, |], || ||] makes it reason- 
able to anticipate the achievement of the predicted BCS- 
transition H Q ^ §]• The existence of a gap in the excita- 
tion spectrum of the superfiuid Fermi-gas will be the first 
issue to address, and several methods for detecting it have 
already been proposed Jl(], [yj. Trapped atomic Fermi- 
gases will allow to study and test fermion-fermion pair- 
ing theories in a tunable, controlled manner. For exam- 
ple the classic problem of the BCS-BEC crossover when 
the interparticle attraction varies (l2) could be studied 
using the possibility of tuning the interatomic scatter- 
ing length using Feshbach resonances 0, |l3). Besides 
the standard superfiuid phenomenology, gases of Cooper- 
paired atoms are expected to have properties which are 
specific to atomic gases only and not present, or not eas- 
ily realizable, in metallic superconductors or Helium. For 
instance, the trapping potential has a major effect on the 
characteristic lengths of the superfiuid Fermi-gas Q . 

In this paper we propose a way to investigate the 
Josephson effect in trapped superfluids of Fermionic 
atoms. We find a phenomenon that is unique to atomic 
Fermi-gases, namely an asymmetry in the Josephson cur- 
rents corresponding to the "up" and "down" spin states. 
We assume that Fermionic atoms in four different hy- 
perfine states (we label them \g), \g'), |e) and |e')) are 
trapped simultaneously in an optical trap — recently all- 
optical trapping and cooling below the degeneracy point 
of the two lowest hyperfine states of 6 Li has been demon- 
strated 0]. The s-wave scattering lengths are assumed 
to be large and negative between atoms in states \g) and 



and 



and the 



\g'), as well as between those in |e) 
chemical potentials [i g ~ fi g ' and /i e ~ /v ■ For all other 
combinations of two atoms in different states the scat- 
tering length is assumed to be small and/or the chemical 
potentials unequal. This leads to the existence of two su- 
perfluids, one consisting of Cooper pairs of atoms in the 
states \g) and \g'), and the other of |e)-|e') pairs. The 



configuration is experimentally challenging, but in prin- 
ciple possible by the choice of right atoms and hyperfine 
states, adjusting the number of atoms, and tuning the 
scattering lengths in magnetic fields by using Feshbach 
resonances J?], [ll| . 

The two superfluids are coupled by driving laser- 
induced transitions between the states \g) and \e) with 
the laser Rabi frequency £1 and detuning 5, and between 
the states \g') and |e') with the Rabi frequency W and de- 
tuning 5'. For a Raman- process these are effective quan- 
tities. If several lasers are used then in order to be able 
to see the Josephson oscillations they should maintain 
their phase coherence for a time much longer than the 
inverse of the detunings. In the case of metals, the two 
superconductors are spatially separated and connected 
by a tunneling junction. In our scheme, the superfluids 
share the same spatial region and are connected by the 
laser-coupling of the atoms' internal states; this resem- 
bles the internal Josephson effect in atomic Bose-Einstein 
condensates jl!) or in superfiuid 3 He-A fll| . 

For metallic superconductors the a.c. Josephson cur- 
rent is driven by applying a voltage over the junction - 
here the role of the voltage is played by the laser detun- 
ings. The difference is that the detunings can be differ- 
ent for the two states forming the pair; in the metallic 
superconductor analogy this would mean having a dif- 
ferent voltage for the spin-up and spin-down electrons, 
a situation which has not been investigated in the con- 
text of metallic superconductors. There is an interest- 
ing connection to recent experiments on superconductor- 
ferromagnet proximity effects, where the chemical poten- 
tials of the spin up and down electrons are slightly differ- 
ent in the ferromagnet due to the exchange interaction 

We consider a system described by the standard BCS- 
theory. The laser interaction is assumed to be a small 
perturbation and its effect is calculated using linear re- 
sponse theory. The observable of interest is the change 



2 



in the number of particles in one of the states, say |e) or 

In the rotating wave approximation the interaction of 
the laser light with the matter fields can be described by 



a time-independent Hamiltonian in which the detunings 
6 and 5' play the role of an externally imposed difference 
in the chemical potential of the two states. The total 
Hamiltonian becomes then H = H + Hj>, where 



Ho = H BCS + L e + 0| dr^ e (f)Mr) + U 9 ~ J dr^ g (f)^ g (r) 
■ + Y) /^,(r)Vv(r) + L g - Yj j dr^ gl (r)^ g ,{r). 
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FIG. 1: The effect of the detunings S and 5' of the two lasers. 
The initial chemical potentials (in the absence of laser cou- 
plings) were fj, g — fi g i and /i e = [i e i . The states (gg') and 
(ee') are Cooper-paired. 



Here /i e and \x g are the chemical potentials of the Fermi 
gases before the laser was turned on, fi g = fi g i and [i e = 
fi e > in order to allow standard BCS pairing. The Hamil- 
tonian Hbcs is the BCS-approximation of the matter 
Hamiltonian with the chemical potential included Jl8| . 
Figure |l| presents a schematic view of what happens in 
this case: the laser detunings shift the chemical poten- 
tials of the four hyperfine states. 

The transfer Hamiltonian is given by 

H T = j dm^ e {r)^ g {r^ + W{^l{r)^ e {r) 

dm'{r)^\,{r)i> g ,{r) + ft'*(f)^,(r)4' (f)(1) 

with fl(r) and fl'(r) characterizing the local strength of 
the matter-field interaction. 

The main observable of interest, the rate of transferred 
atoms from, say, state \g) to state |e), is defined by 

I e = ^ f dr(*(t)\ti(r)Mr)\*(t)) (2) 



(the definition for I e i is similar) and can be further 
evaluated with the help of the Schrddinger equation 



fofa |*(t)> = 8\9(t)) as 



(3) 

In the following we call I e the current in analogy to metal- 
lic superconductors where the flux of electrons out of the 
superconductor constitutes the electrical current. 

We introduce an interaction representation with re- 
spect to H and use linear response theory with respect 
to Ht- Validity of the linear response theory requires 
that the laser intensity is small and the transfer of atoms 
can be treated as a perturbation. 

We split the result for the current I e into a part which 
corresponds to the Josephson current I e j and to the part 
which describes normal single-particle current I e s, I e = 
I e j + I e s- The single-particle current can be evaluated 
at finite temperature using the standard techniques of 
superconducting Green's functions; we present however 
only the result for T = and positive detunings: 



',5' 



= - 2 ^E 



drVt{r)v e n {r)u a m {f) 



5« + 4 n 5). 



Here the triplet (u n ,v n ); e n is a solution of the (nonuni- 
form)_Bogoliubov-de Gennes equations for superconduc- 
tors |l8| and 8 = fi e — fj, g + S. This is the standard Fermi 
Golden rule result and very similar to the ones obtained 
in |pd| . The current I e s is zero when S < 2 A since pair 
breaking is required for single particle excitations. Next 
we concentrate on the Josephson current - this is non- 
zero also for detunings 6 that are smaller than twice the 
gap energy. 

The Josephson current becomes 
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The current J e /j is the same only that <5 and 5' are inter- 
changed. Note that the oscillating term is proportional to 
both of the detunings whereas the rest of the expression 
is proportional only to 5'. For the choice of a homoge- 
neous system (large trap, local density approximation) 
and a constant laser profile the expression simplifies into 



I eJ = I (6')sm[(5 + 6')i\ 
I e ,j = I (6)sin[(5 + 5')i\. 



(5) 
(6) 



Both partners of the pair thus oscillate in phase, with the 
same frequency 6 + 5'. But the amplitudes are different 
whenever the detunings 5 and 5' differ. This means that 
more atoms are transferred, say, in the \g) — |e) oscillation 
than in the \g') — \e') one. 

A simple expression for Iq(5) can be derived when we 
assume identical superfluids, that is A = A' and /i 9 = 
He = M : 



Io(5) 



V2^V A2n2 



-„ VC 2 + A 2 (4C 2 +4A 2 -^ 2 )' 



where V is the volume of the sample and the variable £ 

.2 

is the continuous version of = ^ — fi. Since A <C H, 
the result can also be written as 



Io(5) 



-co N /^ +A 2(4^2 + 4A 2 _ £2) 



A plot of the intensity Iq as a function of the detuning 
5 is shown in Fig. |2|. This result shows a divergence 
at 5 = 2 A, which reflects the divergence of the density 
of states for the two superconductors at the gap. For 
quite a large range of detunings, the amplitude Iq of the 
Josephson current is approximately constant — thus no 
asymmetry effect will be visible. The asymmetry is most 
pronounced when the timescale of the oscillation, that is, 
1/(5 + 5') is close to 1/(2A). Note that 1/(2A) can also 
be understood as the Cooper pair correlation time based 
on the uncertainty principle. 

According to the conventional intuitive picture of the 
Josephson effect, the particles forming a Cooper-pair tun- 
nel "together" through the junction. Therefore our re- 
sult seems counterintuitive at first glance. The physics 
becomes, however, more transparent by a closer look 
at the equation (Q). For simplicity, we consider here 
the transfer process to one direction only, from the su- 
perfluid (gg 1 ) to (ee'), which corresponds to the first 
denominator in (Q). In the initial state \g) is paired 
with \g'), in the final state |e) with |e'). The process 
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FIG. 2: The current 7o (arbitrary units) as a function of de- 
tuning 5, given in units of A. 



has, however, an intermediate state as indicated by the 
second-order form of the observable, and the intermedi- 
ate states corresponding to the observables I e j and I e >j 
are different: For I e j, \g') has been transferred into |e'). 
Therefore its pairing partner \g) is left as an excitation 
in the superfluid (gg 1 ) with the energy e a n and |e') be- 
comes an excitation in the superfluid (ee') with the en- 
ergy t e m = t e m . In contrast, for I e >j, the atom in \g') 
remains as a quasiparticle of the energy e a n = t 9 n in the 
superfluid (gg 1 ) and |e) becomes an excitation in the su- 
perfluid (ee'). For I e j, the initial energy of the Cooper 

pair was — | ) + ( (i g — 4j- J and the energy of the in- 



termediate state is — ■ 

(for explanation, see FigJ^). The relative energy of the 
intermadiate state with respect to the initial state is 

in § 

energy 



e m + 5' , which is precisely the first denominator 
For I e 'j, the initial energy of the pair is again 

<5\ 1 ( ,. <5' S 
2/ 



H g — y ) , but the intermediate state has an 



1 



e^J , or a relative 



- [(^ + |) + 

energy t 9 n + e e m + 5. In summary, the intermediate states 
of the transfer processes for "spin up" and "spin down" 
atoms have different energies and this results in different 
amplitudes for I e j and I e >j. 

The asymmetry in the currents implies the exis- 
tence of excitations in the superfluids. We have ana- 
lyzed the many-body wavefunction of the system in the 
Schrodinger picture and indeed it contains excitations 
corresponding to the asymmetry. Specifically, the analy- 
sis confirms that the so called Fermi surface polarization 
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((N e ) - (N e ,))/((N e ) + (N e ,)) ~ «JV e ) - (N e ,))/({N e ) + 
{N e ')o) is non-zero and oscillates as f(S, 5') cos((<5 + S')t) 
where f(6, S) = 0. We also found that time-independent 
perturbation theory is not sufficient to reveal the asym- 
metry in the amplitudes: the simple treatment of 113] 
applied to our system results in symmetric currents be- 
cause the ansatz used does not allow any excitations. 
This and the fact that the oscillation is most pronounced 
for timescales of the order of the Cooper-pair correlation 
time indicates that the effect is related to the dynamics 
of the superfluid state. 

To observe the Josephson effect one should be able 
to measure the number of particles in two of the states, 
e.g. \e) and |e'), at different stages of the oscillations, 
either destructively or non-destructively. The scale of 
the gap energy is for typical systems l-100kHz, which 
mea ns that the highest time resolution needed should 
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of particles accurately is the more challenging part of 
the observation. In we considered laser probing of 
the superfluid Fermi-gas, where the laser was creating 
excitations in the BCS state. The number of particles 
transferred was directly reflected in the absorption of the 
light. Here one can use similar techniques to detect the 
Josephson oscillations in a simple way. 

In summary, we propose a method to realize Josephson 
oscillations in superfluid atomic Fermi-gases. The cou- 
pling between two superfluids is provided by laser light, 
and the laser detuning plays the same role as voltage over 
metallic superconductor junctions. Detunings that affect 
the two atomic internal states involved in pairing can be 
chosen to be different - this would correspond to different 
voltage for spin-up and spin-down electrons. This leads 
to asymmetry in the oscillation amplitudes of the two 
states. The asymmetry is pronounced when the time- 
scale of the oscillation is the same order of magnitude as 
the Cooper-pair correlation time. This is an effect unique 
to atomic Fermi-gases in the superfluid state. 
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